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Abstract 

The excitation of muhiple surface-plasmon-polariton (SPP) waves of different hnear polarization 
states and phase speeds by a surface-rehef grating formed by a metal and a rugate filter, both of 
finite thickness, was studied theoretically, using rigorous coup led- wave- analysis. The incident plane 
wave can be either p or s polarized. The excitation of SPP waves is indicated by the presence of 
those peaks in the plots of absorbance vs. the incidence angle that are independent of the thickness 
of the rugate filter. The absorbance peaks representing the excitation of s-polarized SPP waves are 
narrower than those representing p-polarized SPP waves. Two incident plane waves propagating in 
different directions may excite the same SPP wave. A line source could excite several SPP waves 
simultaneously. 
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I. INTRODUCTION 



Surface plasmon-polariton (SPP) waves are surface waves guided by a planar interface 
of a metal and a dielectric material. SPP waves find applications for sensing, imaging and 
communication [1, 2]. If the dielectric partnering material is isotropic and homogeneous, 
only one SPP wave — that too, of the p-polarization state — can be guided by the metal- 
dielectric interface at a given frequency [3, 4]. If a periodic nonhomogeneity normal to 
the wave-guiding interface is introduced in the dielectric partnering material, multiple SPP 
waves with different polarization states, phase speeds, and spatial profiles can be guided 
by the metal-dielectric interface. This has recently been shown both theoretically [5-7] and 
experimentally [8-10]. In all of these studies, the dielectric partnering material is also locally 
orthorhombic. 

Very recently, we have solved a canonical boundary- value problem [11] to show that 
multiple SPP waves can be guided even if the dielectric partnering material is isotropic — 
provided that material is also periodically nonhomogeneous normal to the interface. This 
is a very attractive result, because both partnering materials are isotropic and because the 
dielectric partnering material can be fabricated as a rugate filter [12-16]. 

The canonical boundary-value problem does not possess direct practical significance, 
because both partnering materials are assumed to be semi-infinite normal to the planar 
interface. Therefore, we set out to investigate the excitation of multiple SPP waves by 
the periodically corrugated interface of a metal and a rugate filter. This grating-coupled 
configuration [2, pp. 35-41] is popular, when the dielectric partnering material is homoge- 
neous, because it allows the excitation of an SPP wave by a nonspecular Floquet harmonic. 
The interplay of the periodic nonhomogeneity of the dielectric partnering material and a 
periodically corrugated interface is phenomenologically rich [17, 18], and should lead to the 
excitation of multiple SPP waves as different Floquet harmonics. 

The relevant boundary-value problem was formulated using the rigorous coupled-wave 
analysis (RCWA) [19, 20]. In this numerical technique, the constitutive parameters are ex- 
panded in terms of Fourier series with known expansion coefficients, and the electromagnetic 
field phasors are expanded in terms of Floquet harmonics whose coefficients are determined 
by substitution in the frequency- domain Maxwell curl postulates. The accuracy of solution 
is conventionally held to depend only on the number of Floquet harmonics actually used 
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in the computations [21]. The RCWA has been used to solve for scattering by a variety of 
surface- rehef gratings [20-23] , generally with both partnering materials being homogeneous. 

The theoretical formulation of the boundary- value problem is provided in Sec. II and the 
numerical results are discussed in Sec. III. Concluding remarks are presented in Sec. IV. 
An exp(— iwt) time- dependence is implicit, with u denoting the angular frequency. The 
free-space wavenumber, the free-space wavelength, and the intrinsic impedance of free space 
are denoted by = w^/eo/Io, Aq = 27r/A;o, and tjq = a/ /lo/eo, respectively, with jjo and 
eo being the permeability and permittivity of free space. Vectors are in boldface, column 
vectors are in boldface and enclosed within square brackets, and matrixes are underlined 
twice and square-bracketed. The Cartesian unit vectors are identified as u^;, u^, and u^. 
The superscript T denotes the transpose. 

II. BOUNDARY- VALUE PROBLEM 
A. Description 

Let us consider the schematic of the boundary- value problem shown in Fig. 1. The 
regions z < and z > are vacuous, the region < 2; < di is occupied by the dielectric 
partnering material with relative permittivity 6^(2;), and the region d2 < z < hj the 
metallic partnering material with spatially uniform relative permittivity e^. The region 
di < z < d2 contains a surface-relief grating of period L along the x axis. The relative 
permittivity eg(x, z) = eg(x ± L, z) in this region is taken to be as 




(1) 



for 2; G ((ii, ^2), with 





(2) 



and 








1, 



(3) 
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The depth of the surface-rehef grating defined by Eq. (2) is d2 — di. This particular grating 
shape is chosen for the ease of fabrication; however, the theoretical formulation given in the 
remainder of this section is independent of the shape of the surface-relief grating. 




^^ted light 



Surface waves 



Transfr^tted light 
♦ 



FIG. 1: Schematic of the boundary-value problem solved using the RCWA. 



In the vacuous half-space 2 < 0, let a plane wave propagating in the xz plane at an angle 
6 to the z axis, be incident on the structure. Hence, the incident, refiected, and transmitted 
field phasors may be written in terms of Floquet harmonics as follows: 



Hir),c(-'' 

E,e/(r 



neZ 



^70 



r))exp[^(e^x + fc(")^)] ,z<0, 



E (^"^i"^ + Pn4"^) e^P + - 4)] } , ^ > 4 , 

ngZ 

r/o^' {ptt^ - s„4")) exp {z [^^^ + 4"^ - 4)] } , > 4 



(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
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where kx = ko sin 6 + = 27r/L, and 



fci") = <! . (10) 

(k^J^^f - kl kl<{k^^) 



The unit vectors 



and 



(11) 



iL(") /,(") 

= + -7-U. (12) 

/Co Ko 

represent the s- and p-polarization states, respectively. 
B. Coupled ordinary differential equations 

The relative permittivity in the region < z < ^3 can be expanded as a Fourier series 
with respect to x, viz., 

e(x, 2;) = e^*^^ {z) exp{innxx) , z G [0,6/3] , (13) 

where = 27r/L, 



ed{z) , z G [0,4] , 

i/o^^9(^'^)^^' 2;G(rfi,d2), (14) 
Cm , 2; G [(i2, d^l , 



and 



T Cofx, 2;) exp(— mKj;x)(ix , z G [^1,(^2] 
I , otherwise 
The field phasors may be written in the region < 2; < ^3 in terms of Floquet harmonics as 

E(r) = + E'^;\^)^y + E^-\z)iy.] exp(ze^^) 

H(r) = [H^:'\z)v.x + H^;\z)uy + i/i")(^)u.] exp(ze^^) 



^G[0,4], (16) 



n.ez 



with unknown functions ExJ,ziz) and HxJ^ziz). 



5 



Substitution of Eqs. (13) and (16) in the frequency-domain Maxwell curl postulates results 
in a system of four ordinary differential equations and two algebraic equations as follows: 

k^:^^Ej;\z)=koVoHi-\z), 

±Hi!^\z)-tk^:^^Hi-\z) = -^ 
dz r]o 

d_ 

dz 



Y^e^''-'^\z)E^;^\z] 



(17) 

(18) 
(19) 
(20) 



^0 



5^6(— )(Z)E(^ 



^^Hz) 



kt^Hl-\z) 



^0 ^ 



(21) 
(22) 



Equations (17)-(22) hold ^z G (0,^3) and Vn G Z. These equations can be reformulated 
into an infinite system of coupled first-order ordinary differential equations. This system 
can not be implemented on a digital computer. Therefore, we restrict \n\ < Nt and then 
define the column {2Nt + l)-vectors 

[X.(^)] = [Xi-'^^Kz), Xi-^^\zl XW(,), xr-^)(^), Xf')(.)]^ (23) 

for X G {E, H} and cr G {a;, z}. Similarly, we define {2Nt + 1) x {2Nt + l)-matrixes 

[^J = diag[e^] , [g{z)] = [e^--\z)] , (24) 

where diag[/ci"''*] is a diagonal matrix. 

Substitution of Eqs. (19) and (22) into (17), (18), (20) and (21), to eliminate eI"^ and 
gives the matrix ordinary differential equation 
d 



-[i{z)] = z\E{z)]-[{{z)], ze{0,d,), 

where the column vector [i{z)] with 4(2 A^t + 1) components is defined as 



[f(^)]= [E.(^)]", [Ey{z)f, ^o[H.(^)]^ Vo[liy{z)] 



(25) 



(26) 



and the 4(2iVj + 1) x 4(2A^j + l)-matrix [P{z)\ is given by 



M 



P iz) 

[a] 
[a] 

[ol 



(27) 
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Whereas [Q] is the {2Nt + 1) x {2Nt + 1) null matrix and [X] is the {2Nt + 1) x {2Nt + 1) 
identity matrix, the three non-null submatrixes on the right side of Eq. (27) are as follows: 

1 



K 



K 



1 



K 



(28) 
(29) 
(30) 



C. Solution algorithm 



The column vectors [f(0)] and [f((i3)] can be written using Eqs. (4)-(9) as 

Y1 



[f(0)] 



[A] 

m 



=h 



where 



A] 


= H 


-Nt) 


ni-^t+l) 

1 




ni^t-l) iNt) 






-Nt) 


ni-^t+l) 


. . , U-p , 




R] 


= H 


-Nt) 


ri-^t-hl) 

' s ' • 




(7Vt-l) (TVt) 
■■■> ' s ' ' s ' 




rp C 

P 


-Nt) 


'p ) • 


' p ) • 


(7Vt-l) ^(TVt)!'^ 
■■■> ' p ' p J 


[T] 




-Nt) 


A-Nt+l) 


+(0) 


.(7Vt-l) .(TVt) 

■1 ''S 1 ''S 1 




i- 


Nt) 




+(0) 
., , .. 


dNt-l) ANt)!^ 
1 ^p ' P J ' 



and the non-zero entries of {'iNt + 2) x [ANt + 2)-matrixes 



are as follows: 



\ ^ J nm ko 

fy±) = _i 

\ h / nm ' 



n = m + 2A^t + 1 , 
n = m — 2Nt — 1 , 
n = m e [1, 2A/'t + 1] , 
n = me [2Nt + 2, ANt 



2]. 



(31) 



(32) 
(33) 
(34) 

(35) 
(36) 
(37) 
(38) 



In order to devise a stable algorithm [20, 21, 23, 24], the region < z < di is divided 
into A''^ slices and the region di < z < d2 into Ng slices, but the region d2 < z < d^ is kept 
as just one slice. So, there are -\- Ng -\- 1 slices and A^^ + A^^ + 2 interfaces. In the jth 
slice, J e[l,Ns + Ng 

Zj + Zj^i 



-r 1], bounded by the planes z = zj^i and z = zj, we approximate 



P 



z e {zj, Zj_i) 



(39) 
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so that Eq. (25) yields 

[f(z,_i)] = y ^. ■ exp {-zA, J ■ y 7^ • [f (z,)] 



(40) 



where Aj = 2^ — [G] is a square matrix comprising the eigenvectors of as its 

columns, and the diagonal matrix contains the eigenvalues of iii the same order. 

Let us define auxiliary column vectors [T]^. and transmission matrixes . by the relation 
[23] 



J e [0, Nd + Ng + 1] 



where do = 0, 



Na+Ng+l 



Nd+Ng+1 



(41) 



(42) 



To find [T]^. and [£] . for j G [0, Nd + Ng], we substitute Eq. (41) in (40), which results in 



the relation 



(43) 



J e [1, Nd + Ng + l], 

where and \_^- are the upper and lower diagonal submatrixes of [^^.; respectively, 

when the eigenvalues are arranged in decreasing order of the imaginary part. 

Since [T]^. and cannot be determined simultaneously from Eq. (43), let us define [23] 

[n-_i = exp{-A.^;}-^;-[T]^. , 



(44) 



where the square matrix \W\ " and its counterpart \W\ . are defined via 

^3 J 

Substitution of Eq. (44) in (43) results in the relation 

r w 

[exp{-zA,[^]^.}. [a;.{^];}".exp{zA, 



(45) 



(46) 
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From Eqs. (45) and (46), we find in terms of [^]^ _|_^ 

[7] _ L=J 



. After partitioning 



and using Eqs. (31) and (41), [R] and [T]q are found as follows: 



'|T1„" 




'm: - 


Y- 




-1 








. . 




m'o - 















(47) 



[A] 



(48) 



Equation (48) is obtained by enforcing the usual boundary conditions across the plane 
z = 0. After [T]g is known, [T] = [T]j^^_^_j^^_^_-^ is found by reversing the sense of iterations in 
Eq. (44). 



III. NUMERICAL RESULTS AND DISCUSSION 
A. Homogeneous dielectric partnering material 

Let us begin with the dielectric partnering material being homogeneous, i.e., ed{z) is 
independent of z. This case has been numerically illustrated by Homola [2, p. 38] and we 
adopted the same parameters: Aq = 800 nm, = 1.766 (water), = —25 + 1.44z (gold), 
and L = 672 nm. The incident plane wave is p polarized (ctp"^ = 6no 

and ai"^ = OVn G Z) 

and the quantity of importance is the absorbance 

Aft 

A,=i- y: (krf + r + lef + ) R« K^VA^f ] , (49) 

n=-Nt 

which simplifies to 

Nt 

^p=i- E (k?T + i^f ) [ev^f ] , (50) 

n=-Nt 

because all materials are isotropic. 

Figure 2 shows the variation of Ap versus the incidence angle 6 for a sinusoidal surface- 
relief grating defined by [2] 



g[x) = ^(^2 - di) 



1 + sin 



L 



(51) 
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FIG. 2: Absorbance Ap as a function of the incidence angle 9 for a sinusoidal surface-relief grating 
defined by Eq. (51). Black squares represent di = 1500 nm, red circles di = 1000 nm, and blue 
triangles di = 800 nm. The grating depth (^2 — di = 50 nm) and the thickness of the metallic 
layer (^3 — ^2 = 30 nm) are the same for all three cases. The arrow identifies an SPP wave. 
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FIG. 3: Same as Fig. 2 except that the surface-relief grating is defined by Eq. (2) with Li = 0.5L. 



and Fig. 3 shows the same for the surface-relief grating defined by Eq. (2) with Li = 0.5L. 



For computational purposes, we set A^^ 

2 



1, Nn 



the reflectances 



r. 



(n) 



Re 



in) /uiO) 



50, and Nt = 10, after ascertaining that 
2 



(n) 



Re 



ki^^/kf'^ converged 



and the transmittances 
for all n E [-Nt,Nt]. 

Each figure shows plots of Ap vs. 6 for three different values of the thickness di, in order 
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to distinguish [25] between 

(i) surface waves [11], which must be independent of di for sufficiently large values of that 
parameter, and 

(ii) waveguide modes [26, 27], which must depend on di 

as has been shown elsewhere [28, 29]. In both figures, an absorbance peak a.t 9 12° for all 
three values of di indicates the excitation of an SPP wave. 

The relative wavenumbers /ko of a few Floquet harmonics at 6* = 12° are given in 
Table I. The solution of the canonical boundary-value problem (when both partnering 
materials are semi-infinite along the z axis and their interface is planar) [2] shows that the 
relative wavenumber n/ko of the SPP wave that can be guided by the planar gold- water 
interface is 

«/A;o = V^dem/{ed + e™) = 1-3784 + 0.0030^ . (52) 

A comparison of Table I and Eq. (52) confirms that an SPP wave is excited at 6* = 12° as the 
Floquet harmonic of order n = +1. We note that the absorbance peak in Fig. 3 is not only 
wider than in Fig. 2, but also of lower magnitude, which points out the critical importance 
of the shape function g{x) of the surface-relief grating. The incidence angle determined by 
Homola [2, p. 38] is approximately 11°, the small difference between his and our results being 
(i) due to the different methods of computation and (ii) the fact that, while Homola had 
semi-infinite dielectric and metallic partnering materials, we have the two of finite thickness. 

TABLE I: Relative wavenumbers k^'^ /k^ of Floquet harmonics for a gold- water grating when 
6 = 12°. Boldface entries signify SPP waves. 





n = -2 


n = —1 


n = 


n = 1 


n = 2 


e = 12° 


-2.1730 


-0.9826 


0.2079 


1.3984 


2.5889 



B. Periodically nonhomogeneous dielectric partnering material 

Now let us move on to the excitation of multiple SPP waves by a surface-relief grating 
where the dielectric partnering material has a periodic nonhomogeneity normal to the mean 
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plane of the surface-relief grating: 

, z > , (53) 



sm TT- 



where 2f2 is the period. We chose Ua = 1.45 and Ub = 2.32 from an example provided by 
Baumeister [30, Sec. 5.3.3.2 ]. For all calculations reported in the remainder of this paper, 
we chose the metal to be bulk aluminum {emet = ~56 + z21) and the free-space wavelength 
Aq = 633 nm. The surface- relief grating is defined by Eq. (2) with Li = 0.5L. We fixed 
Nt = 8 after ascertaining that the absorbances for Nt = 8 converged to within < 1 % of 
the absorbances calculated with Nt = 9. The grating depth d2 — di = 50 nm and the 
thickness 6^ — 62 = 30 nm were also fixed, as their variations would not qualitatively affect 
the excitation of multiple SPP waves. Numerical results for f2 = Aq and ^2 = 1.5Ao are now 
presented. 



1. n = Xo 

Let us commence with Q = Xq. The solution of the corresponding canonical boundary- 
value problem results in five p-polarized and one s-polarized SPP waves [11], the relative 
wavenmbers n/ko being provided in Table II. To analyze the excitation of s-polarized SPP 
waves in the grating-coupled configuration, we calculated the absorbance 

n=Nt 

As = l- J2 + |^i"^r)Re[A;(")Af] (54) 

n=-Nt 

for ai""* = 6no and ap""* = OVn G Z. Both Ap and Ag were calculated as functions of 6 for 
di G {4Ao, 5Ao, 6A0}, with Ng and A'^^ selected to have slices of thickness 2 nm in the region 
< z < di but 1 nm in the region di < z < d2. Figures 4, 5(a), 5(b), and 5(c) present the 
absorbances as functions of 6* for L = Aq, 0.75Ao, O.TAq, and 0.55Ao, respectively. 

For all three values of (ii, a peak is present at 6' = 37.5° in the plots of Ap vs. 6 in Fig. 4. 
The relative wavenumbers /ko of several Floquet harmonics at this incidence angle are 
given in Table III. At ^ = 37.5°, k^x^/ko = 1.6088 is close to Re[K/ko] = 1.61782, where n/ko 
is the relative wavenumber of an SPP wave in the canonical boundary- value problem [11] 
as provided in Table II; likewise, kx /ko = 1.3912 is close to Re[K/A;o] = 1.36479 (Table II) 
for another solution of the canonical boundary-value problem. Thus, the Ap-peak could 
represent the grating-coupled excitation of either one or two p-polarized SPP waves. To 
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TABLE II: Relative wavenumbers k/Zcq of possible SPP waves obtained by the solution of the 
canonical boundary- value problem [11] for 0, = Xq. Other parameters are provided in the beginning 
of Sec. IIIB. 



s-pol 


1.48639 + 0.00132i 










p-pol 


1.36479 + 0.00169i 


1.61782 + 0.00548i 


1.87437 + 0.00998i 


2.06995 + 0.01526i 


2.21456 + 0.00246f 




9 (deg) 

FIG. 4: Absorbances Ap and Ag as functions of the incidence angle 9 for a surface-relief grating 
defined by Eq. (2) with Li = 0.5L, when Aq = 633 nm, J7 = Aq, and L = Aq- Black squares 
are for di = 6Ao, red circles for di = 5Ao, and blue triangles for di = 4Xq. The grating depth 
(^2 — di =50 nm) and the thickness of the metallic layer (d^ — ^2 = 30 nm) are the same for all 
plots. Each arrow identifies an SPP wave. 

resolve this issue, we computed Ap vs. 6 for L = 1.02Ao (not shown). For the new grating 
period, the Ap-peak shifted to a higher value of 6. This shift showed that the SPP wave at 
the Ap-peak in Fig. 4 is excited as a Floquet harmonic with a positive index n and not the 
one with a negative index. So the Ap-peak represents the excitation of a p-polarized SPP 
wave by a Floquet harmonic of order n = 1 with k^x^ /Uq = 1.6088. 
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TABLE III: Relative wavenumbers /ko of Floquet harmonics when L = Xq. Boldface entries 
signify SPP waves. 





n = -2 


n = —1 


n = 


n = 1 


n = 2 


9 = 28° 


-1.5305 


-0.5305 


0.4695 


1.4695 


2.4695 


9 = 37.5° 


-1.3912 


-0.3912 


0.6088 


1.6088 


2.6088 



A peak is also present at ^ = 28° in the plots of Ag vs. 6, for all three values of di in 
Fig. 4. The solution of the canonical boundary- value problem [11] indicates the excitation 
of an s-polarized SPP wave when k/Zcq = 1.48639 + 0.00132. The real part of this relative 
wavenumber is close to fci^V^o = 1-4695 for 6 = 28°; see Table III. So the Ag-peak represents 
the grating-coupled excitation of an s-polarized SPP wave. 



TABLE IV: Same as Table III except for L = 0.75Ao. 





n = -2 


n = —1 


n = 


n = 1 


n = 2 


9 = 32.5° 


-2.1294 


-0.7960 


0.5373 


1.8760 


3.2040 


9 = 51° 


-1.8895 


-0.5561 


0.7771 


2.1104 


3.4438 


9 = 64° 


-1.7679 


-0.4345 


0.8988 


2.2321 


3.5655 



Since not all possible p-polarized SPP waves (predicted after the solution of the canonical 
boundary-value problem) can be excited with period L = Aq of the surface-relief grating, 
the grating period needs to be changed to excite the remaining SPP waves. The plots of 
Ap vs. 6* for L = 0.75Ao are presented in Fig. 5(a), again for di G {4Ao, 5Ao, 6Ao}. The 
figure shows three Ap-peaks at 6 = 32.5°, 51°, and 64° that are present for all three chosen 
values of di. The relative wavenumbers of several Floquet harmonics at these values of the 
incidence angle are given in Table IV. The Ap-peak at 6 = 32.5° represents the excitation 
of a p-polarized SPP wave as a Floquet harmonic of order n = 1 because kP/ko = 1.8760 IS 
close to Re [1.87437 -|- 0.00998z] in Table II. Remarkably, the Ap-peak at 6 = 51° represents 
the excitation of the same SPP wave as a Floquet harmonic of order n = —2. Finally, the 
Ap-peak at 6 = 64° represents the excitation of another p-polarized SPP wave as a Floquet 
harmonic of order n = 1. 
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FIG. 5: Absorbance Ap as a function of the incidence angle for a surface-relief grating defined by 
Eq. (2) with Li = 0.5L, when Aq = 633 nm, = Aq; and (a) L = 0.75Ao, (b) L = 0.7Ao, and (c) 
L = O.SSAq. Black squares are for di = 6Ao, red circles for di = 5Ao, and blue triangles for di = 4Ao. 
The grating depth (^2 — (ii = 50 nm) and the thickness of the metallic layer (ds — d2 = 30 nm) are 
the same for all plots. Each arrow identifies an SPP wave. 



Two Ap-peaks, at 6 = 27° and 38°, are present for all values of di in Fig. 5(b) for 
L = O.TAq. The relative wavenumbers of several Floquet harmonics at these two values of 
the incidence angle are given in Table V. The Ap-peak at 6* = 27° represents the excitation 
of a p-polarized SPP wave because k^^ /k^ = 1.8826 at this incidence angle is close to 
Re[fi://i;o] = 1.87437 (Table II) of ap-polarized SPP wave predicted by the canonical problem. 
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TABLE V: Same as Table III except for L = O.TAq. 





re = -2 


re = —1 


n = 


re = 1 


n = 2 


e = 27° 


-2.4031 


-0.9746 


0.4534 


1.8826 


3.3111 


e = 38° 


-2.2415 


-0.8129 


0.6157 


2.0442 


3.4728 



The Ap-peak aX 6 = 38° represents another p-polarized SPP wave because kl ' Iko = 2.0442 
is close to Re[K/kQ] = 2.06995 (Table II). 

TABLE VL Same as Table III except for L = 0.55Ao. 





re = -2 


n = —1 


n = 


n = 1 


re = 2 


9 = 25° 


-3.2137 


-1.3956 


0.4226 


2.2408 


4.0590 



All the SPP waves that can be guided by the planar interface of the rugate filter (for 
the chosen value of Q) and aluminum have been shown to be excited by the grating-coupled 
configuration except the one with n/ko = 1.36479 + 0.00169i (Table II). This p-polarized 
SPP wave was found to be excited with a surface-relief grating of period L = O.SSAq. 
The variation of Ap vs. 6 for this case is shown in Fig. 5(c). The Ap-peak at 9 = 25° 
is independent of the thickness di of the rugate filter. The relative wavenumbers k^^ /ko 
of several Floquet harmonics at this value of the incidence angle are given in Table VI. 
Comparing Tables II and VI, we see that ki ^V^o is close to Re [1.36479 + 0.00169i] and 
k^x^/ko is close to Re [2.21456 + 0.00246z]. This suggests that either of the two or both SPP 
waves are excited at 6 = 25°. To resolve this issue, we computed Ap vs. 6 for L = 0.57Ao 
(not shown) and the Ap-peak shifted to a smaller value of the incidence angle. This shift 
indicated that the absorbance peak at 6* = 25° represents an SPP wave due to a Floquet 
harmonic of negative order. So the Ap-peak at 6 = 25° represents the excitation of the SPP 
wave as a Floquet harmonic of order n = —1. 

2. n = l.SAo 

The relative wavenumbers of possible SPP waves that can be guided by the planar in- 
terface of the chosen rugate filter and the metal are given in Table VII for Q = 1.5Ao. 
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In this case, the solution of the canonical boundary-value problem indicated that multiple 
s-polarized SPP waves can also be guided in addition to multiple p-polarized SPP waves. 

Absorbances Ap and A^, calculated for di G {6Ao, 7.5Ao, 9Ao} and 9 G [0°,90°), are 
presented in Figs. 6, 7, 8(a), and 8(b) for L = O.8A0, 0.7Ao, O.QAq, and 0.75Ao, respectively. 
For the computations, the region di < z < d2 was again divided into 1-nm-thick slices; 
however, the region < z < di was divided into 3-nm-thick slices to reduce the computation 
time, after ascertaining that the accuracy of the computed reflectances and transmittances 
had not been adversely affected. 



TABLE VII: Same as Table II except for Q = I.BXq. 



s-pol 


1.61507 + 0.00114 


1.78735 + 0.00078i 




p-pol 


1.40725 + 0.00052i 


1.54121 +0.00374i 


1.71484 + 0.0049i 


1.88541 + 0.00739i 


2.11513 + 0.0045f 


2.02159 + 0.0130H 



TABLE VIII: Same as Table III except for L = O.SXq. 





n = -2 


n = —1 


n = 


n = 1 


n = 2 


e = 17° 


-2.20276 


-0.9576 


0.2924 


1.5424 


2.7924 


e = 21° 


-2.1416 


-0.8916 


0.3584 


1.6084 


2.8584 


e = 27.5° 


-2.0382 


-0.7882 


0.4617 


1.7118 


2.9618 


e = 32° 


-1.9701 


-0.7201 


0.5299 


1.7799 


3.0299 


e = 40° 


-1.8572 


-0.6072 


0.6428 


1.8928 


3.1428 


e = 44° 


-1.8053 


-0.5553 


0.6747 


1.9447 


3.1947 



In the plots of Ap vs. 6 in Fig. 6, the excitation of p-polarized SPP waves is indicated 
at four values of the incidence angle: 9 = 17°, 21°, 27.5°, and 40°. The relative wavenum- 
bers k^x^^ /ko of a few Floquet harmonics at these values of the incidence angle are given in 
Table VIII. The Ap-peak at 6* = 17° represents the excitation of ap-polarized SPP wave, be- 
cause k^x^ /k^ = 1.5424 is close to Re [k/Zcq] = 1.54121 (Table VII), where k/Zcq is a solution of 
the canonical boundary-value problem. The Ap-peak at 6 = 21° also represents ap-polarized 
SPP wave because /i;i~^V^o = 2.1416 is close to Refft/Zco] = 2.11513 (Table VII), which is 
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FIG. 6: Absorbances Ap and Ag as functions of the incidence angle 9 for a surface-relief grating 
defined by Eq. (2) with Li = 0.5L, when Aq = 633 nm, 0, = I.SAq, and L = O.SAq. Black squares 
are for di = 9Ao, red circles for di = 7.5Ao, and blue triangles for di = GAq. The grating depth 
{d2 — di = 50 nm) and the width of the metallic layer (ds — d2 = 30 nm) are the same for all the 
plots. Each arrow indicates an SPP wave. 

another solution of the canonical problem. The relative wavenumber /ci^V^o = 1-6084 at 
6 = 21° was ruled out, after examining the plots for L = 0.82Ao (not shown). Similarly, the 
peak at ^ = 27.5° is due to the excitation of another p-polarized SPP wave as a Floquet 
harmonic of order n = 1 and not of the order n = —2. At 6 = 40°, both k^x^ /ko = 1.8928 
and ki~^^/ko = 1.8572 are close to Re [1.88541 + 0.00739?] (Table VII). However, the plots 
of Ap vs. 6 for L = 0.82Ao (not shown) made us conclude that the SPP wave is excited as 
the Floquet harmonic of order n = 1. 

In the plots of As vs. 6 in Fig. 6, two peaks at ^ = 32° and 44° are present for all 
values of di. The Ag-peak at 6 = 32° represents the excitation of an s-polarized SPP wave 
because fci^V^o = 1-7799 is close to Re [k/Z^o] = 1-78735 (Table VII), which is a solution of 
the canonical boundary-value problem for an s-polarized SPP wave. The other A^-peak at 
9 = 44° represents the excitation of the same s-polarized SPP wave as a Floquet harmonic 



18 



of order n = —2. 




FIG. 7: Same as Fig. 6 except for L = 0.7 Xq. 



TABLE IX: Same as Table III except for L = O.IXq. 





n = -2 


n = —1 


n = 


n = 1 


n = 2 


9 = 10° 


-2.6835 


-1.2549 


0.1736 


1.6022 


3.0308 


9 = 16° 


-2.5815 


-1.1529 


0.2756 


1.7042 


3.1328 


9 = 27.5° 


-2.3954 


-0.9668 


0.4617 


1.8903 


3.3189 



For L = 0.7Ao, the absorbances Ap and Ag are presented as functions of 6 in Fig. 7. The 
relative wavenumbers k^"^ /ko of several Floquet harmonics at those values of the incidence 
angles where peaks are present independent of the value of di are given in Table IX. In the 
plots of Ap vs. 6, the peak at 6* = 16° represents the excitation of a p-polarized SPP wave 
because k^x^ /ko = 1.7042 is close to Re[K,/ko] = 1.71484 (Table VII). Similarly, the Ap-peak 
at 6 = 27.5° represents the excitation of a p-polarized SPP wave, as k^^^/ko = 1.8903 is close 
to Re[K/ko] = 1.88541 (Table VII). In the plots of A^ vs. 9, the peak at 6* = 10° represents 
the excitation of an s-polarized SPP wave because k^^ /ko = 1.6022 at this angle is close to 
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FIG. 8: Absorbance Ap as a function of the incidence angle 9 for a surface-relief grating defined 
by Eq. (2) with Li = 0.5L, when Aq = 633 nm, Q = I.SAq; and (a) L = 0.9Xq and (b) L = 0.55Ao- 
Black squares are for di = 9Ao, red circles for di = 7.5Xq, and blue triangles for di = 6Ao. The 
grating depth {d2 — di = 50 nm) and the width of the metallic layer (da — ^2 = 30 nm) are the 
same for all the plots. Each arrow indicates an SPP wave. 



TABLE X: Same as Table III except for L = O.OAq. 





n = -2 


n = —1 


n = 


n = 1 


n = 2 


e = 43° 


-1.5402 


-0.4291 


0.6820 


1.7931 


2.904 


e = 51.5° 


-1.4396 


-0.3285 


0.7826 


1.8937 


3.0048 


e = 55° 


-1.4031 


-0.2920 


0.8192 


1.9303 


3.0414 



In Fig. 8(a), three peaks are present for all values of di at 6 = 43°, 51.5°, and 55°. The 
relative wavenumbers of several Floquet harmonics at these values of the incidence angle are 
given in Table X. The comparison of Tables VII and X shows that each Ap-peak represents 
the excitation of a p-polarized SPP wave. Similarly, each of the three Ap-peaks at 6* = 22.5°, 
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TABLE XI: Same as Table III except for L = 0.75Ao. 





n = -2 


n = —1 


n = 


n = 1 


n = 2 


e = 22.5° 


-2.2840 


-0.9506 


0.3827 


1.7160 


3.0414 


e = 33.5° 


-2.1147 


-0.7814 


0.5519 


1.8885 


3.2186 


e = 46° 


-1.9473 


-0.6140 


0.7193 


2.0526 


3.3860 



33.5°, and 46° in Fig. 8(b) represents the excitation of a p-polarized SPP wave as a Floquet 
harmonic of order n = 1. The relative wavenumbers of a few Floquet harmonics at these 
values of the incidence angle are given in Table XI. At each angle, fci^V^o is close to the real 
part of one of the k/Zcq of the solutions of the canonical boundary- value problem problem 
provided in Table VII. 

C. General conclusions 

In the last two subsections, we have deciphered a host of numerical results and identified 
those absorbance peaks that indicate the excitation of SPP waves in the grating-coupled con- 
figuration, when the dielectric partnering material is periodically nonhomogeneous normal 
to the mean plane of the surface-relief grating. We found that 

(i) the periodic nonhomogeneity of the dielectric partnering material enables the excita- 
tion of multiple SPP waves of both p- and s-polarization states; 

(ii) fewer s-polarized SPP waves are excited than p-polarized SPP waves; 

(iii) for a given period of the surface-relief grating, it is possible for two plane waves with 
different angles of incidence to excite the same SPP wave (Figs. 5(a) and 6); 

(iv) not all SPP waves predicted by the solution of the canonical problem may be excited 
in the grating-coupled configuration for a given period; 

(v) the absorbance peaks representing the excitation of p-polarized SPP waves are gener- 
ally wider than those representing s-polarized SPP waves; and 

(vi) the absorbance peak is narrower for an SPP wave of higher phase speed (i.e. smaller 
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Let us note that some other combination of the periodic functions ^di^) and g{x) may allow 
all solutions of the canonical boundary-value problem to be excited in the grating-coupled 
configuration. 

The solution of the canonical boundary- value problem [11] indicates that the period 2Q 
of the rugate filter needs to be greater than a certain value in order for more than one SPP 
waves to be excited, and the excitation of s-polarized SPP waves to exist requires an even 
larger period. However, the number of possible SPP waves increases as the period increases 
up to a certain value. We chose Q = Xq and Q = 1.5Ao to allow the excitation of multiple 
p-polarized SPP waves for both values, and multiple s-polarized SPP waves for the second 
value. Our numerical results confirm that the conclusions on the number of SPP waves 
drawn in the predecessor paper [11] also hold for the grating-coupled configuration. 

IV. CONCLUDING REMARKS 

The excitation of multiple surface-plasmon-polariton (SPP) waves by a surface-relief grat- 
ing formed by a metal and a dielectric material, both of finite thickness, was studied theoret- 
ically using the rigorous coupled-wave-analysis technique for the practically implementable 
setup. The presence of an SPP wave was inferred by a peak in the plots of absorbance vs. 
the angle of incidence 6, provided that the 6'-location of the peak tuned out to be indepen- 
dent of the thickness of the partnering dielectric material. If that material is homogeneous, 
only one p-polarized SPP wave, that too of p-polarization state, is excited. However, the pe- 
riodic nonhomogeneity of the partnering dielectric material normal to the mean plane of the 
surface-relief grating results in the excitation of multiple SPP waves of different polarization 
states and phase speeds. In general, the absorbance peak is narrower for an s-polarized SPP 
wave than for of a p-polarized SPP wave, and the absorbance peak is narrower for an SPP 
wave of higher phase speed. 

Since the electromagnetic field radiated by a line source can be considered as a spectrum 
of plane waves propagating at all angles [31, Sec. 2.2], the grating-coupled configuration 
discussed in this paper can be used to excite multiple SPP waves simultaneously by a line 
source. The excitation of multiple SPP waves may be significant for practical applications — 
for example, to increase the absorption of light in solar cells due to the increased possibility 
of excitation of SPP waves [32]. This application is currently under investigation by the 



22 



authors. 
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